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Abstract 

We have used coarse-grained simulation methods to investigate the effect of stretching-induced 
structure orientation on the proton conductivity of Nation-like polyelectrolyte membranes. Recent 
experimental data on the morphology of ionomers describe Nafion as an aggregation of polymeric 
backbone chains forming elongated objects embedded in a continuous ionic medium. Uniaxial 
stretching of a recast Nafion film causes a preferential orientation of these objects in the direction 
of stretching. Our simulations of humid Nafion show that this has a strong effect on the proton 
conductivity, which is enhanced along the stretching direction, while the conductivity perpendicular 
to the stretched polymer backbone is reduced. Stretching also causes the perfluorinated side chains 
to orient perpendicular to the stretching axis. This in turn affects the distribution of water at low 
water contents. The water forms a continuous network with narrow bridges between small water 
clusters absorbed in head-group multiplets. 
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I. INTRODUCTION 



In their role as proton-conducting membranes, ionomers are an important component of 
many hydrogen fuel cells. In these materials, the interplay among the short-range inter- 
actions between the hydrophobic backbone polymer and the hydrophilic terminal groups, 
and the long-range Coulomb interactions between the electrostatic charges on the terminal 
groups and the protons induces a nanophase separation into proton-rich and proton-poor 
domains. A general model for the phase morphology of ionomers has been proposed by 
Eisenberg et al.—, according to which a few head groups combine to form multiplets that 
restrict the mobility of the backbone chain segments directly attached to them. A spher- 
ical geometry is often assumed for these multiplets, whose sizes are typically less than a 
nanometer 2 . The average distance between these multiplets is mostly dictated by the con- 
centration of head groups relative to that of backbone monomers. These multiplets then 
form microdomains, and one observes a microphase separation that can serve to facilitate 
proton diffusion. Recent experimental data on the morphology of ionomers describes the 
aggregations as elongated objects embedded in a continuous ionic medium 3 . A goal of much 
ionomer research is to increase the proton conductivity, and hence to make membranes that 
can operate under very low humidity conditions, as a higher conductivity results in a higher 
output power density in fuel cells containing these membranes. 

Proton conduction itself is a complex process, which strongly depends on the thermal and 
mechanical history of the membrane. The mechanical history involves the manufacture of 
the membrane, which is usually achieved by one of two common procedures: solution casting 
or extrusion. The former method is used to make membranes from a solution of dissolved 
ionomer by allowing the solvent to evaporate from the solution. While this technique is 
suitable only for small-scale laboratory production, it has the advantage producing isotropic 
membranes with no residual preferential orientation of their backbone within the plane of 
the membrane. 

Most commercially available ionomeric membranes are fabricated by extrusion of a 
molded sample. This leads to a preferred orientation of the ionomer backbone 4 , the extent 
of which depends on the draw speed of the extrusion. This inherent structural anisotropy of 
the backbone matrix is believed to be a reason for the susceptibility to tearing or cracking 
of the membrane in any swelling or drying processes. This creates technical problems in 
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keeping the membrane taut in the changing temperature-humidity conditions encountered 
in fuel cells. An equally important issue concerns the effect of this anisotropy on proton con- 
ductivity. While it has generally been assumed that such effects are significant, the problem 
of predicting the consequences of mechanical strain remains largely unsolved. 

There have been several experimental studies in which membranes were uniaxially 
stretched in order to probe the effects of strain on internal morphology. Gebel et al.— ana- 
lyzed the form of ionic domains in unstretched and stretched membranes, and showed that 
mechanical stretching induces ordering in the ionomer backbones. Elliot et al^- detected 
the anisotropy in stretched membranes by performing scattering experiments. Barbi et al.— 
used membrane elongation measurements to confirm the classical ionomer domain model, 
in which inverted micelles are interconnected by channels. As expected, uniaxial stretching 
of recast Nation causes a preferential orientation of the Nation backbone in the direction 
of stretching, and this is morphologically similar to the anisotropy in extruded membranes. 
Cable et al.—, for example, stretched Nation and found the proton conductivity to be higher 
in the plane of the membrane than normal to it. A slightly different result was found by 
Lin et al.—, who noted little change in transverse conductivity on stretching, but found an 
improved fuel cell performance, as compared to Nation 117 and unstretched recast Nation, 
as a consequence of a lowered methanol permeation rate. In a related experiment, Elabd 
et al— measured the conductivities of an ionic block copolymer, and observed appreciable 
anosotropy. Oren et al.— found that the conductivities of membranes could be rendered 
anisotropic through the alignment of suspended particles, with higher conductivity in the 
alignment direction. These studies reveal the significant impact that organized and oriented 
structures can have in increasing proton conductivity. While there have been several inten- 
sive simulational studies of the effects of naturally occurring anistropies on the transport 
properties of Nafion-like membranes^^yL^JiyLii, there are no studies of which we are aware 
of the effects of mechanical stretching on an equilibrated sample. 

In this paper we investigate numerically to what extent induced orientation in the back- 
bone polymer affects proton diffusion and conductivity in Nafion-like ionomers. The nature 
of our united-atom model for the ionomer is described in Section II, and the details of our 
simulation procedure are given in Section III. Section IV reports our results for stretch- 
induced morphological changes in the ionomer, while results for the proton conductivity 
are given in Section V. In Section VI we discuss our conclusion that stretching of solvated 
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membranes leads to anisotropy of the conductivity, with the stretch-induced ordering of the 
backbone increasing the longitudinal protonic conductivity while decreasing the transverse 
conductivity. 



II. MODEL 

Because it is necessary to model quite large volumes of material in order to see the types of 
morphological changes in which we are interested, it is unfortunately not feasible to perform 
fully atomistic simulations. We thus employ the same united-atom model for Nafion^^^ 
as was used in our previous work—. Within this approach the CF2 and CF3 groups of the 
backbone and sidechains, the ether oxygens of the sidechain, the sulfur atom and the O3 
oxygen group of the sulfonates are modeled as Lennard- Jones (LJ) monomers of diameter 
a = 3.5A. All the ether oxygens and fluorocarbon groups are assigned zero partial charges. 
The electrostatic charges are located entirely on the SO3 sulphonate head groups and on 
the H + protons. 

The total potential energy of the system is then given as 



Here the first two terms on the right-hand side are pair interactions between nonbonded 
monomers. The first term is the Lennard- Jones interaction 



where Elj = 0.2 kcal/mol, and rjj = |r< — fj\ is the separation distance between monomers 
% and j. At a temperature of 300K, this value of Elj would correspond to an energy of 
about ksT/3. The parameter a is 1 for hydrophobic-hydrophobic (HH) interactions, and 
0.5 for both hydrophobic-hydrophilic (HP) and hydrophilic-hydrophilic (PP) interactions 
between monomers. In the PP case, the LJ potential was modified to be purely repulsive by 
truncating it at its minimum, where r = 1.122c, and raising it by the addition of an amount 
Elj- These interaction parameters were chosen to agree in most instances with the Nafion 
model of Paddison et al.—. 

The electrostatic interaction between charged sulphonate groups and protons is 



^total — f^LJ + U q + £/bond + Wangle + ^dihedral- 



(1) 




(2) 



(3) 



i>j 
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Here q\ (I = i,j) is +1 for protons, +1.1 for sulfur atoms, and —2.1 for the combined oxygens 
of the sulfonate groups, thus ensuring that the total charge of a sulfonate head group SO3 is 
— e. The distance- dependent dielectric permittivity e(r) in Eq. ([3]) reflects the fact that, at 
small separations, there is no intervening material between any two charges. There is then 
no screening of the Coulomb interaction, and the effective permittivity is close to unity. At 
modest distances the charges will be separated by matter that may be water, ionic groups, 
or backbone polymer, and the energy of interaction will be reduced by the shielding effect of 
this material. The interaction between sulphonates, both inside the sulphonate head group 
clusters surrounded by polymeric matrix, and between two different sulphonate clusters 
separated by polymeric matrix, is thus a complicated function that includes the effects of 
image charge distributions at the boundaries^. When water is present, dielectric saturation 
and water immobilization effects in the hydrophilic phases of the ionomeric membrane 2 ^ 2 - 
add complications. While the best approximation for e(r) is still a matter for debate 2 ^, 
it is clear that the permittivity must increase with distance. We choose the form e(r) = 
1 + e#((r I ' — l)/(r/<7 + l)) 10 , which smoothly increases from e = 1 near a chosen charge 
to being the dielectric permittivity of the bulk at large distances. Varying the value of 
€b between 8 and 20 has only a negligible effect on all the observed nanophase separations 
in the system under consideration. The effective dielectric permittivity e e g, calculated as a 
volume average of e(r) over a typical multiplet domain, appears to be in good agreement 
with reported dielectric constants for ionomers in Ref.— . 

The last three terms on the right side of Eq.flTJ represent the potential energy of the 
bonded segments of the molecules. The two-body bond-stretching potential 

U hond (R) = \ E h{R-Ro) 2 (4) 

all bonds 

has a simple Hookean form with stretching constant fc& = 700 kcal/mol A 2 (roughly 1200 
/cgT/A 2 at room temperature) and unstretched bond length i? =l-54 A (0.44a). For the 
three-body angle-bending potential, 

W«) = | E k e (e-e )\ (5) 

all bond pairs 

we use an equilibrium bending angle #0 = 110° and a bending force constant kg = 120 
kcal/mol deg 2 (or about 200 /c^T/deg 2 ). Finally, the dihedral (four-body) component of 
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the total energy is written as 

^dihedral {4>) = ^ Yj ^ (1 - d COs(30) ) , (6) 

all bond triplets 

with the parameters having values d — — 1 (+1) and fc^=18.1 kcal/mol ~ 10.8 kgT (6.2 
kcal/mol ~ 3.7fceT) for backbone (sidechain) segments, respectively. 

We consider two models for the ionomer. In the first model, which is a highly-simplified 
representation of a completely dry membrane, each proton is permanently attached to a 
sulfonate group. Each side-chain head-group is thus electrostatically neutral, and consists 
of two dipoles, connected head to tail. We refer to this as the bound-proton model. 

In the second, and more realistic, model, the proton is more loosely bound to its sulfonate 
counter-ion by a Coulomb attraction. We refer to this as the solvated model. The activa- 
tion energy needed to overcome the electrostatic attraction between the oppositely charged 
H + and SO3 ions, and thus detach the proton, depends on the average distance between 
sulfonates and on the water content A in the membrane. This parameter A is defined as 
the number of water molecules per sulfonate head group, and varies between and 25 for 
solvated membranes. It is believed that A > 5 is sufficient to loosen the bonds tying the 
protons to the head groups and to achieve full dissociation. These free protons contribute to 
the conductivity principally through diffusion and through the Grotthus mechanism. When 
A < 5, hopping diffusion, in which protons make a transition from one sulfonate to an 
adjacent one, becomes the main contributor to the conductivity. 

We use a simple point charge (SPC) fluicr^ 1 ^ to model the water molecules in our sim- 
ulations. These move in a medium having a distance-dependent dielectric permittivity e(r) 
that we assume to be independent of A. This includes the case of A = 0, where in the 
absence of SPC water molecules the proton conductivity occurs only through the hopping 
mechanism. 

III. SIMULATION DETAILS 

The simulations were performed in stages. The first step was to grow a continuous 
backbone with attached pendant sidechains by using Monte-Carlo techniques. The bond 
length, and the bending and dihedral angles for this initial polymer were fixed to be Rq, 6q 
and 0o respectively. After molecular-dynamics (MD) runs of a few picoseconds, which were 
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used to equilibrate the system, as this usually needed to respond to the strong steric repulsion 
of overlapping monomers, the constraints on angles and bonds were removed. Following this, 
the side chains were detached from the backbone skeleton in an approach previously adopted 
by Vishnyakov 2 ^, and the backbone skeleton was cut into segments 2 ^. This greatly increases 
the relaxation rate of the ionomer, and allows a rapid equilibration of the system. The system 
then consists of N s = 1000 sidechain segments and iV& = 1400 backbone segments. Taking 
into account that the sulfonic acid groups are hydrophilic while the ethers and fluorocarbon 
groups are hydrophobic^, we adopted the following coarse-grained representation. The side 
chain architecture is written as 7H+3P for a dry membrane and 7H+2P for a humidified 
membrane. Here the letters H and P stand for hydrophobic and hydrophilic monomers 
correspondingly. The backbone segments are fully hydrophobic with a 14H architecture. 

Molecular-dynamics runs with a Langevin thermostat were performed for the segmented 
polymer for time periods up to 50ps in a constant-iV^T ensemble. We imposed periodic 
boundary conditions for a cube of side L = 30er and used the Lekner summation method 2 ^ 
to handle the long-range electrostatic interactions between charged particles. In order to 
verify that the system had not become trapped into a metastable glassy state, we repeated 
each run with several different initial configurations. 

In the next stage of the simulations, the segments were reassembled back into the branched 
chain characterizing the original Nation. This was achieved by the simultaneous introduc- 
tion of the bonds and angular constraints described by Eqs. dl|-(j6]) between the ends of 
each backbone segment, which united them into a single chain. Similar bond and angular 
constraints between the tail monomer of each detached side chain and the median monomer 
of every backbone segment connected the side chains to the backbone polymer. To avoid 
the formation of star-like branched polymers, only a single occupancy of the backbone at- 
tachment sites was permitted. The unequal numbers of backbone segments and side chain 
segments resulted in the formation of a polymer in which there were varying numbers of 
backbone monomers between the points of attachment of the side chains. The simulation 
was then resumed, and run until a new equilibrium was achieved. Once the system was fully 
equilibrated, the statistically averaged quantities of interest were gathered during the next 
2-3ns of run time for the initial case of unstretched membranes. 

The stretched samples were constructed by the application of a set of forces designed to 
emulate the effects of a longitudinal stress. The continuous backbone polymer was considered 
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as a set of Nb connected segments, and a weak force f± was applied to the end monomer of 
each segment. The direction of this force was in the positive z direction for the end with the 
greater value of its z coordinate, and in the negative z direction for the end with the lesser 
value of its z coordinate in an extended- zone scheme where the periodic boundary conditions 
were not applied. This force tends to align the backbone segments in the z direction without 
having significant influence on the side-chain dynamics. The strength of f± was restricted to 
values much smaller than the forces arising from the bond, angle, and dihedral potentials. 

It usually took 50ps to equilibrate the stretched membrane. Following this, runs of 
2-3ns were performed to gather statistically significant results for the morphology of the 
stretched material. For studies of the proton conductivity, a weak electric field was applied 
to the equilibrated stretched membrane, and additional 50ps runs were performed to reach 
steady-state conditions before Non-Equilibrium MD (NEMD) simulation runs of 2-3ns were 
undertaken. In this way the effects of both stretching and applied fields could be determined. 

IV. RESULTS FOR STRETCHED-MEMBRANE MORPHOLOGY 

A snapshot of an equilibrated membrane with a water content of A = 3 molecules per 
head group and in the absence of an applied stretching force is given in Fig. [T](a). The 
hydrophobic backbone and the side chains are shown as lines (red and blue respectively in 
the online version). The spheres indicate the sulfonate head groups of the sidechains, and 
show cluster-like aggregations of different sizes and shapes, in agreement with experimental 
observations 29 . A moderately stretched membrane with dimensionless stretching force f=l 
is shown in Fig. [T^b). The polymer backbone shows ordering in the direction of the applied 
force, which is parallel to the z axis. We determine the induced backbone ordering Sb and 
sidechain ordering S s by calculating the order parameter 



an approach commonly used in liquid crystal systems. Here j = b,s for the backbone 
polymer and sidechains respectively, the angular brackets (...) denoting statistical averaging 
over the NVT ensemble. The orientational angle 9 for each segment is defined as the angle 
between the vector connecting the first and last monomers of the segment and the applied 
force f , as illustrated in Fig. [2j The backbone matrix is considered as a set of Nb segments 
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of 14 monomers each. If 9 is larger than tt/2, the order of monomer numbering along the 
segment is reversed for a proper accounting its orientation. In a similar manner we define the 
orientational angle 9 for sidechains as also shown in Fig. [2J Here 9 is the angle between the 
stretching force f and the vector pointing from the first monomer of the side chain towards 
its terminal group. 

The order parameters Sb and S s , calculated from Eq. ([7]), are shown in Fig. [3j As ex- 
pected, membrane stretching leads to ordering in the ionomer conformation. The backbone 
ordering, shown in Fig. [3]^a), is greatest in the bound-proton model and decreases in the 
solvated case as the water content A increases for a fixed force f . Figure [3](b) reveals the com- 
pletely different response of the side chains to the stretching force. When the bound-proton 
membrane is stretched, the sidechains tend to orient themselves parallel to the stretching 
force. However when a solvated membrane is stretched, the sidechains are more likely to 
be perpendicular to the stretching force, in accord with previous experimental results^. 
A possible explanation for these different responses to stress lies in the effect of the hu- 
midity on multiplet formation. In the absence of water, the dipolar interactions between 
sulfonates leads to the formation of dense multiplets containing 10-20 terminal groups, as 
seen in Fig. H](a). The Coulomb forces of attraction between these dipoles limit the amount 
of molecular relaxation of the sidechains in response to stretching. The backbone defor- 
mation is thus followed by the rearrangement of multiplets via elongation in one direction 
and shrinking in the other direction. This is evident from a comparison of the curves for 
stretched and unstretched membranes in Fig. Hh- 

In contrast, when the membrane is solvated, the water content lowers the densities of 
the multiplets. These then coalesce to form larger clusters in which a mixture of solvent 
molecules and conducting protons gives the head groups more dynamic flexibility. As a 
result , the sidechains can achieve a gain in entropy by rearranging themselves perpendicular 
to the direction of the applied stress. This is the reason why only a slight deformation of the 
multiplet geometry for the solvated membrane is seen in Fig. H](b). In summary, the order 
parameter of the membrane depends on its hydration level, the absorbed solvent being an 
obstacle for the backbone ordering Sb, but facilitating the side-chain ordering S s . 

The most important aspect of the effect of stress on morphology concerns the formation 
of channels through which proton transport might occur. An indication of this is shown 
in Fig. [51 in which the water distribution in a stretched ionomer is displayed. The water 
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molecules congregate in elongated clusters in association with the distribution of side-chain 
head-groups, and are oriented in the direction of the applied stress. During the course of the 
simulation, the clusters create temporary bridging connections allowing the water to diffuse 
easily back and forth along the stretching direction. This morphology and dynamics of the 
water undoubtedly facilitates proton diffusion parallel to the stretching direction. 

In order for proton conduction to occur, it is necessary for there to be a continuous 
network of pathways through the hydrophilic material in a membrane. This makes it useful 
to have a quantitative measure of the connectivity of the hydrophilic channels. To this end 
we have developed a method^ aimed at detecting hydrophilic pathways in a network of 
hydrophilic clusters. The only input parameter for this method is the maximum separation 
distance parameter r m between neighboring particles in the pathway. Our method calculates 
the distribution function P{Ls) of the hydrophilic channel length L$, which is defined as 
the distance between the initial and end points of the pathway. We applied this method to 
determine how the stretched-membrane morphology alters the channeling in the hydrophilic 
network of clusters. The distribution function P(Ls) when r m is chosen to be 2a (which 
is 7 A, the average separation between sulfonates in Nation) is shown in Fig. El both for 
stretched and unstretched membranes. Both curves have a similar short-range behavior and 
show a maximum at a most probable channel length equal to about 3a. This is roughly the 
average size of the sulfonate multiplets in solvated membranes, as was seen in Fig. H^b). The 
long-range behavior of P(Ls), however, shows a dependence on the mechanical treatment of 
the membrane. In a stretched membrane the hydrophilic-sulfonic channels in the membrane 
are 50% longer than in an unstretched membrane. We note that the calculated P(L S ) 
is an instantaneous measure of the pathways in a network of clusters. Because of the 
constantly changing patterns in the distribution, which are caused by the diffusion of ions 
and water and the motion of head groups, the real pathways are effectively longer, and can 
even traverse the whole system. Consequently one expects that the protons in a stretched 
membrane can continuously travel along the stretching axis z. We explore the resulting 
protonic conductivities in the next section. 
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V. RESULTS FOR THE CONDUCTIVITY OF A STRETCHED MEMBRANE 



In this section we analyze the effect of membrane stretching on the protonic current. 
The electrical conductivity x due to proton current in a solvated membrane can be ob- 
tained by performing NEMD simulations to find the current j in the presence of an external 
electrostatic field E through the relation 

j = x(E)E. (8) 

Here the induced current density j is given by 

N s 

j = (iA0E^,e (9) 
i=i 

where i^e is the E-field component of the velocity of the i th ion, and Ns is the total number 
of protons in the system. 

Because a strong applied electric field may itself have an effect on the morphology, we 
isolate the effects of stretching by examining only the limit Xo °f as E — > 0. The 

longitudinal and transverse conductivities were evaluated by applying fields parallel and 
perpendicular to the stretching direction, and extrapolating the resultant conductivities to 
find their values at zero applied field. The calculated values of Xo f° r the case where E 
is along the membrane stretching axis are shown in Fig. [7] as a function of the membrane 
order parameter <S& for different solvation parameters A. It is seen that ordering in the 
membrane backbone increases the conductivity, with a dependence of Xo(Sb) on 5*6 that is 
close to being linear for all water contents A. Conversely, when we examine the transverse 
protonic conductivity Xo across the stretched membrane, which is shown in Fig. [BJ we find 
that the conductivity decreases with increasing backbone ordering. These two trends can 
be understood as a direct consequence of the channel-like nanophase separation along the 
stretching direction. While the protons can freely travel along these channels as a conse- 
quence of the water cluster bridgings, there are no continuous pathways for them to travel 
across the stretching direction. This anisotropy in the conductivity is most pronounced in 
the dryer samples. 
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VI. DISCUSSION AND CONCLUSION 



The simulation results presented in the previous sections show that the uniaxial stretching 
of a Nafion-like ionomer should be expected to have a significant effect on its proton conduc- 
tivity. The conductivity along the stretching direction is considerably higher than its value 
in an unstretched membrane, while the proton diffusion perpendicular to the stretching axis 
is strongly reduced. When a water-containing sample is stretched, the perfluorinated side 
chains exhibit a tendency to be oriented perpendicular to the stretching force direction. The 
sulphonate clusters deform as the stretching is increased, and show a long-range ordering 
in their distribution across the membrane. This affects the distribution of absorbed water, 
which forms a continuous network with narrow bridges. 

The effect of the water is different for the backbone and sidechain ordering parameters. 
As seen from Fig. [3j water is an obstacle to backbone ordering Sf, parallel to the direction 
of elongation. However, water enhances the sidechain ordering S s perpendicular to the 
stretching direction. For a given stretching force f , the order parameter S& of the membrane 
backbone depends on the solvation parameter A. 

The simulations reported here were carried out at room temperature, and under condi- 
tions of constant stress, and we made no attempt to investigate how the anisotropy in the 
protonic conductivity depends on the temperature at which the membrane was stretched. 
Experimental measurements of crystallinity, however, indicate that membrane stretching at 
elevated temperatures can lead to different results from ambient-temperature stretching^ 1 ^. 
It appears that uniaxial extension of membranes at temperatures above the boiling tem- 
perature of water results in a material with an oriented morphology that persists after the 
samples are removed from the extensional stress. Also, higher-temperature compression- 
molded ionomers exhibit a backbone crystallization when stretched^. It would thus be of 
interest to pursue further studies to enlarge on those presented here. 

Finally, we note that with regard to the issue of improving fuel-cell performance by using 
stretched membranes, our results indicate that the observed increased proton conductivity 
occurs only in the direction of stretching. As a practical matter, the desired direction 
of proton transport is across the membrane, and so the production of a membrane with 
enhanced conductivity due to stress may be difficult to arrange. 
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Figure 1: (Color online) A typical snapshot of the simulation box for (a) unstretched (/=0) and 
(b) stretched (/=1) Nation membrane with water content A = 3. The stretching force / in (b) is 
directed along the z axis. Colored beads represent the end-group oxygens of side chains. Pendant 
side chains and neutral backbone polymer are drawn by lines (blue and red, respectively, in online 
version). Different bead colors correspond to different bead altitudes, with blue for beads at the 
bottom of the simulation box and red for beads at the top of the box. The size of all structural 
elements is schematic rather than space filling. 
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Figure 2: (Color online) This schematic illustration explains the calculation of order parameters 
Sb for backbones (left picture), and S s for sidechains (right picture). The continuous black line 
represents the backbone polymer of the ionomer. Shorter gray lines (green in online version) 
represent the sidechains. The dashed lines are vectors constructed to show the orientation of 
backbone segments and sidechains, with 8 being the angle between the orientation vector and 
stretching direction. 
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Figure 3: (Color online) Backbone order parameter Sb (a) and sidechain order parameter S s (b) as 
functions of applied stretching force / for the bound-proton model and for different water contents 
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Figure 4: (Color online) This figure shows the effect of strain on the sulfur-sulfur pair distribution 
function g ss (r) for (a) the bound-proton model and (b) the solvated ionomer with A = 5. Dashed 
arrows show the average size of multiplets. 
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Figure 5: (Color online) Instantaneous distribution of water clusters in a stretched, solvated mem- 
brane for stretching force / = 1 and water content A = 5 water molecules per sulfonate. 
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Figure 6: (Color online) The distribution function P(L S ) of sulfonate channel lengths L s for un- 
stretched (/=0) and stretched (/ = 1) solvated membranes having a water content A = 5. 
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Figure 7: (Color online) Protonic conductivity xo of a stretched membrane along the stretching 
direction as a function of the backbone order parameter Sb for different membrane solvations. The 
water-content parameters A are indicated below the corresponding lines. 
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Figure 8: (Color online) Protonic conductivity xo of a stretched membrane perpendicular to the 
stretching direction as a function of the backbone order parameter Sb for different membrane 
solvations. The water-content parameters A are indicated below the corresponding lines. 
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